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Abstract 

A detailed description of an original method used to measure the luminosity 
accumulated by the HERA-B experiment for a data sample taken during the 2002- 
2003 HERA running period is reported. We show that, with this method, a total 
luminosity measurement can be achieved with a typical precision, including overall 
systematic uncertainties, at a level of 5% or better. We also report evidence for the 
detection of (5-rays generated in the target and comment on the possible use of such 
delta rays to measure luminosity. 



1 Introduction 



D A precise determination of tlie luminosity is required for the measurement 
of absolute cross sections. The integrated luminosity (£) is defined by 

Np 

>C=— (1) 

CTp 

where Np is the number of events of a given process and dp is the correspond- 
ing cross section. In the case of HERA-B, which is a forward spectrometer [1,2] 
experiment, operated at the 920 GeV proton beam of the HERA accelerator 
at the DESY Laboratory in Hamburg, the proton beam is bunched and in- 
teracts with a nuclear target placed on the halo of the beam. The number 
of proton-nucleus (pA) interactions per bunch crossing is subject to statisti- 
cal fluctuations. For HERA-B, as for all other experiments having a bunched 
beam, the luminosity can be expressed as: 

C = ^^SLlA (2) 

a 

where A is the average number of interactions per bunch crossing BX, Nbx is 
the number of beam bunches crossing the apparatus and a is the interaction 
cross section (for a more detailed discussion, see Sections 4 and 5). As a 
consequence, given the cross section of proton-nucleus interactions, the lumi- 
nosity can be measured by determining A and Nbx- The average number of 
interactions per BX can be obtained from a fully unbiased sample of events in 
various ways: by looking at inclusive quantities which are proportional to the 
number of interactions in one event (such as the number of tracks or the en- 
ergy released in a calorimeter), by counting the number of primary vertices or 
by counting the number of empty events. The first method has the advantage 
of entailing only a rather straightforward analysis of the data, but the signal 
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corresponding to a single interaction must be evaluated precisely and detector 
stability becomes a relatively critical issue. In the second method, the vertex 
reconstruction efficiencies must be known precisely as well as the probability 
of erroneously merging or splitting primary vertices during reconstruction. In 
the third method, the distribution of the number of interactions per bunch 
crossing must be cither known or assumed and the efficiency for detecting 
non-empty events and the impact of noise events must be evaluated. 
After careful studies the HERA-B Collaboration has decided to exploit the 
method based on counting events with evidence of at least one interaction 
(which is equivalent to the third method listed above), since this method 
minimizes the systematic error on the luminosity determination allowing to 
achieve a final precision of about 5%. 

The paper is organized in the following way. In Sections 2 and 3 the main 
features of the HERA accelerator relevant for this analysis and the HERA- 
B detector are briefly described. Section 4 summarizes all of the published 
proton-nucleus cross section measurements which are used for the luminosity 
determination. In Section 5 and 6, the relevant relations for the determination 
of the luminosity are described. In Section 7 we discuss the systematic uncer- 
tainties and comment on delta ray production, while in Section 8 we report 
the results obtained for the interaction trigger (deflned below) data sample. 



2 The HERA accelerator and the target 

HERA is a double storage ring designed for coUiding a 920 GeV proton beam 
with a 26 GeV electron beam. Four interaction regions exist: two of them 

house the general purpose ep detectors HI and ZEUS, while the other two 
accommodate the fixed target experiments HERA-B and HERMES. In the 
following we describe the beam parameters and the filling scheme used during 
the HERA-B data taking period 2002-2003. 

The typical proton current is 80 mA, distributed over 180 bunches with a 
typical bunch length of 1-2 ns. The proton bunches are organized into 3x6 
trains of 10 consecutive bunches each, separated by one empty RF bucket. The 
detailed flUing scheme is shown in Fig. 1. In total there are 220 RF buckets 
with a spacing of 96 ns including a gap of 15 empty buckets at the end to 
provide for a secure beam dump. The average rate of filled bunch crossings is 
8.52 MHz. 

The target system [3] consists of two stations of four wires each. The wires are 
positioned above, below, and on either side of the beam and are made from 
various materials including carbon, titanium and tungsten. Both titanium 
and tungsten targets are wires with a diameter of 50 iim, whereas the carbon 
target is a flat ribbon, 100 /im perpendicular and 500 /im along the proton 
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Fig. 1. Schematic representation of the bunch structure of a HERA proton-ring fill. 

beam. The stations are separated by 40 mm along the beam direction. The 
wires are positioned individually in the halo of the stored proton beam and 
the interaction rate for each inserted wire is adjusted independently. Any 
number of wires can be operated simultaneously. The luminosity measurement 
described herein applies exclusively to single wire runs. 

The steering of the target wires requires a fast and reliable system to provide a 
counting rate proportional to the interaction rate up to the highest interaction 
rates envisaged in the HERA-B design (40 MHz). This is achieved by limiting 
the acceptance of the scintillation counters used to detect interactions to ~ 
10~^. Stepping motors with a nominal step-size of 50 nm controlled by a 10 Hz 
steering loop provide a stable interaction rate. 

One additional complication is that a fraction (typically a few percent) of 
interactions not correlated to any bunch [4] was present. These interactions 
are due to so-called coasting beam protons which have left the separatrix, 
but are still circulating inside the machine, forming a component of the beam 
halo. Based on test measurements, the coasting beam can be regarded as 
a DC-current. The fraction of coasting beam depends on the position of the 
target and the history of the individual proton fill, thus requiring an individual 
correction for each run. As described in Section 7.6, the relevant information 
can be derived from events triggered by a pseudo-random generator. 



3 The HERA-B detector and the data sample 

The HERA-B experiment is a forward magnetic spectrometer with an accep- 
tance extending from 15 to 220 mrad horizontally and to 160 mrad vertically. 
This large angular coverage allows studies in kinematic regions not accessible 
to previous fixed-target high energy experiments. A top view of the detector 
is shown in Fig. 2. The first part of the spectrometer is devoted to tracking 



6 




Fig. 2. Top view of the HERA-B detector. 



and vertex measurements and consists of the target, a silicon vertex detector, 
a magnet and a tracking system. The second part is focused on particle identi- 
fication and includes a Ring Imaging Cherenkov detector, an electromagnetic 
calorimeter and a muon detector. 

The vertex detector (VDS) [5] is placed between the target and the magnet 
and divided in 8 stations. Each station consists of four "quadrants" equipped 
with two double-sided silicon microstrip detectors (50 x 70 mm^, 50 /im pitch) 
each. This system provides a primary vertex resolution of ~ 500 /im along 
the beam direction and ax,y ~ 50 /im in the transverse plane. 
A dipole magnet with a 2.13 Tm field-integral is positioned before the 
main tracking system. Each tracking station consists of several planes of 
MSGC/GEM chambers placed near the beam pipe (Inner Tracker, ITR) [6] 
and several planes of Honeycomb Drift chambers which cover the rest of the 
acceptance (Outer Tracker, OTR) [7]. The detector segmentation is designed 
to cope with the particle flux variation with the distance from the beam pipe. 
Typical momentum resolutions of 5p/p ~ 1% are achieved. 
The particle identification of charged tracks (protons, kaons, etc.) is provided 
by a Cherenkov detector (RICH) installed downstream of the magnet. A /3 ~ 1 
particle traversing the RICH detector produces an average of about 33 hits [8]. 
The electromagnetic calorimeter (ECAL) [9], which provides electron pretrig- 
ger seeds and e/vr separation, is installed after the RICH and the tracking 
system. The ECAL is a Shashlik sampling calorimeter with Pb or W as ab- 
sorber and scintillator as active material. In order to follow the steep radial 
dependence of the particle density, the calorimeter has been structured in three 
sections (Inner, Middle and Outer) with differing granularities. The Muon de- 
tector (MUON) [10] provides the muon pretrigger seeds and the muon identifi- 
cation, and is located in the most downstream part of the detector. It consists 
of four superlayers embedded in an iron loaded concrete absorber. The sen- 
sitive area close to the beam pipe is covered by pixel chambers, while in the 
rest of the acceptance, tube chambers are used. 

The flexibility of the trigger system [11] allows the implementation of a large 
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variety of trigger configurations. The methods described in this paper have 
been used to determine the integrated himinosity of the interaction trigger 
(lA) data sample. The lA trigger selects events with at least one inelastic 
interaction in the target, by requiring either that the RICH has more than a 
minimum number of hits (20) or that the ECAL has more than a minimum 
(1 GeV) energy deposition. The total collected statistics is about 220 million 
events, with an average data acquisition (DAQ) rate larger than 1000 Hz. 
During the data acquisition a sample of randomly trigger events (Zero-Bias) 
was acquired in parallel to the lA trigger at a rate of few Hz, allowing the pos- 
sibility to check the trigger acceptance and stability. Moreover the Zero-Bias 
event sample has been extensively used in the luminosity determination. Due 
to the fact that the same data stream was used both for the various physics 
analyses and for the determination of the recorded luminosity, the dead time 
of the DAQ system cancels exactly and can thus be ignored. 



4 The cross sections 

The total pA cross section atot can be divided into elastic (uei) and inelastic 
{(Tinei) contributions: 

Ctot — CTel + CTinel — CTel + Cmb + CTtsd + CTbsd + dd ■ (3) 

In this context, the cross section a^i is regarded as the sum of the elastic 
{pA — > pA) and quasielastic contribution [pA — > pA*). The inelastic cross 
section includes a minimum bias part {mh) and a diffractive part which can 
be further subdivided into target single diffractive {tsd, pA — > pY), beam 
single diffractive {bsd, pA — > XA) and double diffractive {dd, pA — > XY) 
contributions. 

The values for the total and inelastic cross sections reported in Table 1 were 
obtained using the method of [12] with one exception: in order to minimize a 
possible systematic bias, we do not apply the A" scaling law adopted there. 
Instead we use the experimental results on carbon and tungsten nuclei for 
the total cross section of [13] and for the inelastic cross sections of [14,15,16]. 
All of these measurements were obtained at beam momenta ranging from 
180 to 400 GeV/c and have to be scaled to 920 GeV/c with the prescription 
given in [17]. Due to the absence of data on titanium, the A"' scaling law 
is applied only to interpolate the Al and Fe data of the experiments quoted. 
The elastic cross sections are obtained using Equation 3. The single diffractive 
cross sections are taken from [12]. The experimental results can be compared 
to a theoretical calculation performed in the framework of the Glauber-Gribov 
theory [18]. Both total and inelastic cross sections agree well within 5%, while 
the diffractive contributions exhibit larger discrepancies. As suggested by [19], 
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an average of Model III and Model IV of [18] is used for this comparison. 
The double diffractive part, being neglected in [12], is taken from [18]. 
The quoted errors cover the difference between both models. The minimum 
bias cross section is derived by subtracting all diffractive contributions from 
the inelastic cross section with an error given by the quadratic sum of the 
component errors. 



cross section (mb) 
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W 




351.6 ± 4.0 


1045.± 30. 


2913. ± 43. 


^inel 
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1788. zb 23. 


CTel 


100.9 ± 4.8 


362.5± 30.5 


1125. ± 49. 


Cbsd 


8.4 ± 1.7 


21. 7± 4.3 


41.2 ± 8.2 




9.2 ± 2.3 


14.9± 3.8 


23.9 ± 6.3 


(^dd 


0.7 ± 0.3 


l.lib 0.5 


1.5 ± 0.6 


<^mb 


232.4 ± 3.9 


644.8± 8.1 


1721. ± 26. 



Table 1 

Summary of the proton nucleus cross sections for C, Ti and W nuclei. Details are 
given in the text. 

The total trigger efficiency can be expressed as 

^el ■ CTel + ^mb ' 0"mfe + ^bsd ' '^bsd + ^tsd ' (^tsd + ^dd ' dd , 

Slot = , (4) 

O'tot 

where Ex is the efficiency for triggering on process x. The trigger efficiencies 
are determined from Monte Carlo simulation. FRITIOF 7.02 [20] is used to 
generate minimum bias events in pA interactions, while diffractive events are 
generated by PYTHIA 5.7 [21] which, however, has the disadvantage of not 
taking into account nuclear effects. Nonetheless, since diffractive cross sections 
are small compared to the minimum bias cross section, they contribute little 
to the recorded sample and nuclear effects for diffractive events can be safely 
neglected. 

The detector response is simulated by the GEANT 3.21 package [22]. Realistic 
detector efficiencies, readout noise and dead channels are taken into account. 
The simulated events are processed by the same reconstruction codes as the 
data. The resulting interaction trigger efficiencies are summarized in Table 2. 
It has been checked that the results for the pp processes do not depend on the 

^ Both models are based on the saturated form of the dipole cross section and 
provide a more realistic description compared to Models I and II which assume a 
quadratic dependence. 
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wire position. The small increase of the minimum bias efficiency with increas- 
ing atomic mass number is correlated to the increasing track multiplicity. 



Process 
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£tsd 


PYTHIA 


0.370 
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PYTHIA 


0.578 








^mb 


PYTHIA 


0.941 










FRITIOF 




0.933 


0.953 


0.970 


^tot 


Eq. 4 




0.642 


0.607 


0.586 


Ka 


Eq. 5 




0.960 


0.969 


0.978 



Table 2 

lA trigger efficiencies for the various processes. Ka is defined in Equation 5. 

From these numbers we can conclude that the elastic contribution is negligi- 
ble and that the diffractive processes are suppressed. The dominance of the 
minimum bias part can be illustrated by calculating its detectable fraction Ka 

j-^ Cmfe ■ ^mb 

Ka = (oj 

<^tot ■ ^tot 

given in Table 2. The impact of the uncertainties on luminosities and trigger 
efficiencies will be discussed in Section 7.2. 



5 General remarks on the luminosity determination 



In the following, the luminosity given by Equation 2 will be expressed in terms 
of the total number of events satisfying the lA trigger (Nja), the average 
number of interactions per bunch crossing (Xtot), the trigger efficiency per 
single interaction (etot) and the total hadronic cross section (atot)- In order to 
do this, two assumptions are made: 

• the number of interactions per filled bunch can be described by a single 
Poisson distribution P{n, Xtot), for all bunch crossings in a given data run: 

\n p-\tot 

P{n, Xtot) ^ , , (6) 
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• the trigger efficiency for n interactions, {etot)n, is given by: 

ietot)n = 1 - (1 - etotr (7) 

where Stot is the trigger efficiency for single interaction. 

A test of the vahdity of the first assumption is discussed in Section 7.4, while 
the second assumption has been checked in Monte Carlo studies and by check- 
ing the dependence of the measured Xtot with interaction rate as measured by 
the target steering scintillator hodoscopes (see Section 7.5). 

With these assumptions, the total number of recorded triggers resulting from 
interactions in the target, Nja, is given by: 

oo 

NiA ^Nbx-T. (^(^' ^tot) • = Nbx ■ (l - e-^*--^*-) (8) 

n=0 

where Nbx is the total number of BXs considered. From this equation, given 
the general relationship of Equation 2, we finally obtain 

Because the product £ • A is typically ^ 10% for our data taking conditions, 
the measured luminosity is, to first order, inversely proportional to the trigger 
efficiency and the cross section, while the average number of interactions, Xtot, 
enters only as a second order correction. 

In Equation 9, Nja can be expressed as a function of the number of recorded 
triggers {Ntape) and of the number of background events {N^kg): 

NiA = Ntape - Nbkg = Ntape " (1 " fbkg) (10) 

where fbkg is the fraction of background events in the sample (see Section 7.6). 

Since, as discussed in Section 4, the recorded event sample is dominated by 
minimum bias interactions, the luminosity can be expressed as a function of 
minimum bias quantities 

Ntape ' (1 fbkg) ' Xmb /1 1 \ 

-^io* = ~T\ -E u\ .^ ^^^> 

(1 — e ^m6-^m6j . (7^^ 

where Ka is defined in Section 4 and A^;, is defined in Section 6. This is the 
final expression which will be used to determine the luminosity for each run. 
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6 The determination of X^b 



The determination of Xmb relies on the pseudo random trigger data sample 
acquired in parallel to the I A trigger and the large-acceptance detectors which 
constitute the spectrometer. Specifically, Xmb is obtained by combining the 
information from a variety of subdetectors to also provide a cross-check of 
the stability of the result and the systematic uncertainties due to the detector 
response and the event model of the Monte Carlo. Only the filled bunches of 
HERA were considered. 

The average number of I A per BX can be evaluated with respect to any 
subdetector observable X which depends hnearly on the interaction rate, by 
exploiting the following definition: 



where Nx is the number of events with observable X above a certain threshold 
and e{X) is the corresponding efficiency (i.e. probability that an interaction 
will result in X being above threshold) as evaluated from the FRITIOF simu- 
lation. To avoid possible confusion, we note that X^b as defined by Equation 12 
is close to but not equal to the average number of minimum bias interactions 
per BX. With this definition. Equation 11 is nonetheless exact. The sample 
of random trigger events for all runs is sufficiently large that the statistical 
error is always negligible compared to the systematic error estimate. The list 
of the observables X used for the determination of X^b is given in Table 3. 
As can be seen, two subdetectors are directly involved in this method, namely 
RICH and ECAL, while the VDS and the OTR are indirectly involved when 
the number of reconstructed tracks is considered. 

It is important to note that for the determination of Xmb, no reconstructed 
quantity associated only with the VDS is used. The reason for this is that all 
such quantities were found to be sensitive to the presence of 5-rays generated 
by the proton beam in the target, as will be further discussed in Section 7.7. 

The A values calculated with Equation 12 as a function of two of the X ob- 
servables (namely hrich and e{ECAL, inner)) defined in Table 3 arc shown 
in Fig. 3. These plots show a common feature: the existence of a broad sta- 
tionary point for A^;,. The left plot in Figure 3 shows the Xmb values obtained 
as a function of the threshold applied on the number of hits (hrich) seen in 
the RICH detector. The rise at small values of RICH hits is due to noise 
in the detector, while the smooth increase for large number of RICH hits is 
mainly due to the fact that the Monte Carlo does not precisely reproduce the 
RICH hit multiplicity per event, although the resulting value of A^;, is nearly 
independent of the threshold over a wide range. This trend is confirmed, in a 
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total energy deposition in full ECAL (GeV) 


nclus(ECAL) 


number of reconstructed electromagnetic clusters in ECAL 


hecal 


number of hit towers in ECAL 


ntra 


number of reconstructed tracks (VDS+OTR) 



Table 3 

Description of the eight X quantities used to determine Xmb- 
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Fig. 3. Two typical distributions for the measured average number of interactions 
per bunch crossing as a function of the cut on different reconstructed quantities for 
the tungsten target wire. The distributions for RICH (number of reconstructed hits 
per event , -hrich-), ECAL (total energy released per event in the inner section of 
the calorimeter, -e{ECAL, inner)-) are shown. 

more or less pronounced way (see e.g. the right plot of the same Figure), also 
for the other variables hsted in Table 3. 

As a consequence, for each observable, as the best estimate of X^b the value 
{^min(yX)) is taken to be its minimum value. The best evaluation of Xmb is 
then defined as: 

Xmb = 2^ ^ U3j 

X=l,8 ° 

The values of Xmin{^) obtained from all are in good agreement and their 
spread is used as a measure of the systematic uncertainty. 
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7 Systematic uncertainties and checks 

7.1 General Considerations 

According to Equation 11 and the assumption on which Equation 6 is based, 
the following systematic uncertainties must be taken into account: 

• the uncertainty on K^. arising from the Monte Carlo {MC ) event model 
and the poorly known observation probability of diffractive processes; 

• the uncertainty on the determination of Xmb', 

• the uncertainty associated with deviations of the interaction probability 
distribution from the assumed Poisson distribution (e.g. due to the uneven 
filling of bunches) ; 

• the uncertainty associated with possible out-of-time interactions (i.e. coast- 
ing beam interactions) and fake triggers from detector noise. 

No systematic uncertainty due to reconstruction efficiency appears on this 
list since possible systematic biases due to reconstruction are included in the 
systematic uncertainty assigned to Xmh- 

The determination of Ctot could, however, be biased if the on-line trigger docs 
not operate according to expectations. Possible triggering errors are checked 
for by using special bits written into the event record to indicate the trigger 
decision. The online trigger requirements are imposed offline on the sample of 
random trigger events taken with each I A trigger run and compared to the 
the online decision. It is found that there is no significant inefficiency from the 
online trigger while the percentage of spurious triggers due to electronic mis- 
behavior is typically at level of a few per thousand. Some additional masking 
of noisy channels is found to be necessary, although the effect on the trigger 
acceptance is neghgible. We conclude that the trigger performed according to 
expectations and introduces no additional biases on the measurement of Ctot- 

The general expression for the squared relative uncertainty on the luminosity 
follows from Equation 11: 

/5£tot y _ 5^ ^ ^ Q^(^_^ (14) 

V Ctot ) Ka Xmb \^ Jbkg \^ J Pais ' 

which is the quadratic sum of the relative systematic uncertainties on Ka (see 
Section 7.2), on A^^ (see Section 7.3), on the background and on the Poisson 

assumption (see Equation 6). The last two sources of systematic uncertainty 
will be discussed in Sections 7.4 and 7.6 respectively. Finally, to separate the 
measurement uncertainties from the uncertainties on the present knowledge of 
the total cross section {Ka term), the following quantity (which will be used 
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in Section 8) is defined: 



=X^© -TT ®(-) . (15) 
V -^tot / cjei Arnb \ ^ / bkg \ ^ / Pais 



7.2 Uncertainty on Ka 



The relative uncertainty (7^) on the detectable fraction Ka (defined in Equa- 
tion 5) depends on the uncertainties of the cross sections, quoted in Table 1, 
and on the trigger efficiencies of the various production processes. Given the 
poor knowledge of the structure of final states produced by diffractive pro- 
cesses, the trigger efficiencies are assumed to be fully unknown but limited to 
the range from to 1. Thus, an error of cte^^^ = a^^^^ — = is assigned. 
The error on Smb is not included here as discussed in Section 7.1. The resulting 
uncertainties are summarized in TaV)le 1. 





Ka 


SKa 
Ka 


c 


0.960 


0.023 


Ti 


0.969 


0.018 


W 


0.978 


0.016 



Table 4 

Detectable fractions and their relative uncertainty for carbon, titanium and tungsten. 



7.3 Uncertainty on Xmb 



The method used to determine Xmb is influenced by the Monte Carlo descrip- 
tion of the HERA-B detector as well as the event model of the event generator. 
The resulting uncertainty on Xmb is taken to be the rms spread of the X^b 
values calculated with Equation 13. The typical values obtained are 

r 

^-0.04, (16) 

Xmb 

or better, depending on the target material. 

Possible sources of systematic uncertainty related to the way Xmb is determined 
have been investigated. For example, examination of the observables given in 
Table 3 shows that six of them involve ECAL and are thus possibly subject 
to correlated systematic effects while RICH and VDS-I-OTR appear only with 
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one variable each. For this reason an alternative quantity A^^ is defined as: 

= . (IT) 



The relative discrepancy of \mb and X'^^^ is given by: 




\ — \' 

method 



Xmb 



The mean value of the distribution of this quantity is found to be statistically 
compatible with zero and its rms width is ~ 0.006 or better, depending on 
target material. The smallness of this term, compared to the overall systematic 
uncertainty on A„i5 (see Equation 16), shows that the calculation of X^i, is 
insensitive to the relative weights given to the various methods. 



7.4 Uncertainty on the distribution of the number of I A per BX 



One important assumption is that the number of interactions per HERA ma- 
chine bunch follows a Poisson distribution (see Equation 6). In order to eval- 
uate the systematic uncertainty associated to the non-Poisson behavior of the 
BX population, the total luminosity can alternatively be measured for each 
run, also as a sum over all the BX contributions, i.e.: 



(^mb i=i ( 1 — e ^"ib^^^bij 

The quantity {Ctot)Bx then compared with the total luminosity calculated 
according to the basic procedure (see Equation 11). In this way, we can define 
the systematic uncertainty due to the non Poisson behavior of the beam as: 



_ ^tot {^tot)BX ^20) 

Pais 

The mean of the distribution of this quantity for all the runs with more than 
3 • 10^ events shows a slight shift (~ 0.4%) toward negative values. The rms 
of the distribution is ~ 0.009 or better, depending on target material. 
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1.5 The dependence of X^b on the target interaction rate 

The dependence of Xmb, as determined by the method of Section 6, on the 
target steering hodoscopc rate {Rhod) was checked. In general the dependence 
should be linear at sufficiently low interaction rate. 

For each of the three target materials (C, Ti and W), X^b was determined for 
hodoscope interaction rates of 0.3, 0.5, 1, 3, 5 and 10 MHz. For each target, 
the resulting Xmb values, excluding the 10 MHz point, were fit to a straight 
line. The results of the fit for carbon and tungsten target wires are shown in 
Fig. 4. The 10 MHz point lies below the fit line indicating possible saturation 
of the hodoscopes or possibly a breakdown of the Poisson assumption at high 
interaction rates. The normalized of the linear fit is about one or better 
(for all three set of runs) and indicates a linear relation between interaction 
rates measured by two very different techniques for rates up to 5 MHz. This 
in turn supports the two assumptions made in Section 5 since the hodoscope 
rates do not rely on the Poisson assumption and compensating non-linearities 
in the two methods are unlikely. 




2.5 5 7.5 10 2.5 5 7.5 10 

RhoXMHz) R.od(MHz) 



Fig. 4. Test of linear dependence of the Xmb values determined with the method 
described in Section 6 with respect to the target interaction rate measured by the 
scintillator hodoscope system for the carbon (left plot) and tungsten (right plot) 
target wires. 



7.6 Background estimate 

As shown in Equation 10, the number of recorded events must be corrected for 
background events either from fake triggers (e.g. electronic noise in the RICH) 
or from coasting beam interactions (beam-gas interactions and background 
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from interactions upstream of the target arc negligible). The best method to 
determine the fraction of background events is to compare empty and filled 
bunches using the random trigger events. For this purpose the interaction 
trigger requirements are applied offline to the random trigger events which 
are equally distributed over all 220 bunches. Normalizing to the 180 filled 
bunches, we obtain 

w _ 180 ( N^apejempty BX) \ 

^"'^ 40 \Ntape(filledBX)J " ^ ^ 

A small (approximately 5%) correction is made to account for "in-time" coast- 
ing beam interactions which are well synchronized in time with the detector's 
integration gates and thus are no different from ordinary interactions from 
bunched beam protons. 

The resulting values for fh^g from the random trigger sample are summarized 
in Table 5 and used in Equation 11. To estimate the uncertainty on f^kg, 
the software trigger thresholds are varied over a wide range. The uncertainty, 
SCbkg: is obtained by dividing the difference of the extreme values by ^/T2. The 
uncertainties due to the in-time coasting beam correction discussed above are 
negligible compared to the uncertainties given in the Table. 





fbkg 


(SC\ 
\ C / bkg 


c 


0.031 


0.018 


Ti 


0.057 


0.023 


W 


0.026 


0.019 



Table 5 

The fraction of background events fbkg for each target wire and the relative system- 
atic uncertainty. 

7.7 Production of S -rays in the target 

The presence of (5-rays in the data sample (see Section 6) is both a nuisance, 
since it compromises the VDS based methods, and an opportunity for a sys- 
tematic check of the luminosity calculation, since the luminosity can be esti- 
mated from the observed rate of 5-ray production. The results of a study of 
5-ray production applied to a run taken with the carbon target sample are 
presented and compared to the luminosity estimates given in Section 6. With 
further development, the techniques presented here could be used for a precise 
luminosity determination in experiments using thin targets. 
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The luminosity for a fixed-target experiment in a proton beam is proportional 
to the sum of target path lengths of all protons (A^tot) which traverse the 
target: 



A 



Ntot 



A 



■ E (22) 

where A is the atomic mass of the target material, Na is Avogadro's number, 
p is the target density in (g/cm^) and Zi is the length of the target traversed 
by the ith proton. 

The number of 5-rays (Ns^prod) produced in a kinetic energy (T) interval from 
Tmin to Tmax IS proportional to the same summed target length [17]: 

Ns,prod = 0.154-p ^ ■ / 0-154^;^ E ^^ (23) 

where Z, A, and p are the atomic number, atomic mass and density (in g/cm^) 
of the target and T is in MeV. Tj^ax is approximately 475 GeV for 920 GeV 
incident protons. 

Combining Equation 22 and Equation 23 results in the following equation re- 
lating luminosity to the number of produced 5-rays: 



^ ^ATmin -^A^min ^5,obs fOA^ 

~ 0T54Z ■ ''^'■"'^ ~ 0l54Z ■ ^ ^ 

where Ng^obs is the number of observed 5-rays and £5 is the average probability 
that the S-ray escapes the target and is reconstructed. 

The (5-ray detection efficiency (ss) is evaluated by Monte Carlo. The 5-rays 
are generated according to Equations 27.5 and 27.6 of [17] with a minimum 
kinetic energy threshold of 1 MeV and tracked through the target and de- 
tector using the GEANT3-bascd [22] HERA-B simulation program. Since 
the VDS has acceptance for tracks from the target in the polar angular 
interval 0.01 ^0 ^ 0.7rad, corresponding to a 5-ray momentum range of 
1-88 ^P5^ 10220 MeV/c (see Equation 27.6 of [17]), the 1 MeV kinetic energy 
threshold corresponds to (5-rays which are well outside the detector acceptance. 
The generated Monte Carlo events are subjected to the same reconstruction 
and analysis code (see below) used for the data. 

An average efficiency of ^ 7% after all cuts is found. We estimate a 15% 
relative systematic uncertainty on this number coming from uncertainties in 

the material distribution in the vertex detector and from sensitivity to Monte 
Carlo parameters, in particular to the minimum kinetic energy cutoff for track- 
ing by GEANT, nominally set to 30 keV. (These sources of systematic error 
could in principle be greatly reduced by a more precise inventory of detector 
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materials and by a more thorough study of the tracking of very low momentum 
electrons.) 

Candidate (5-rays are reconstructed using the standard HERA-B VDS recon- 
struction software applied to a pseudo-random-triggered carbon target run and 
then searched for in events from filled bunch crossings which do not pass the 
I A trigger condition. Distributions of track parameters and derived quantities 
from the 5-ray Monte Carlo and data are found to be in close agreement when 
segments which extrapolate to near the average vertex position of hadronic in- 
teractions are removed. Such tracks are typically high momentum tracks from 
hadronic pN interactions. 

We define the impact parameter of a track as the difference between the aver- 
age position of vertices from hadronic interactions (and therefore the average 
impact point of the beam on the wire target) and the track's position when 
extrapolated to the Z-position of the target. Figure 5(a) shows the X-view 
impact parameter {Xip) distribution of reconstructed VDS segments in non- 
lA events which contain a single reconstructed segment with \Xip\ > 1 mm 
which originates in the first VDS layer. The cut on Xip removes a signal from 
high-momentum tracks from hadronic pN interactions (approximately 1/4 of 
the removed tracks form a narrow peak above the relatively broad distribution 
shown in the figure). The data, indicated by the histogram, and the Monte 
Carlo, indicated by the points with error bars, are in good agreement. Note 
that the width of the Xip distribution is largely determined by multiple scatter- 
ing in the VDS and therefore depends on the momentum of the reconstructed 
tracks. The close match between data and Monte Carlo implies that the mo- 
mentum spectra of reconstructed tracks in data and the 5-ray Monte Carlo 
are similar. Figure 5(b), showing the distribution of the X-view impact points 
of the same tracks used in Figure 5(a) at the first VDS layer, also illustrates 
the good agreement between data and Monte Carlo. The corresponding distri- 
butions in the Y-view also agree well with each other as do the distributions 
of track polar angles. The overall close agreement strongly suggests that the 
observed tracks are indeed caused by 5-rays originating in the target. Further 
evidence that the observed tracks are associated with beam protons traversing 
the target comes from the greatly reduced rate of such tracks in empty bunch 
crossings: ~ 3% of the rate in filled bunch crossings. 

The observed rate of 5-ray candidates is 0.068 ±0.003 per BX. When the same 
event and track selection criteria are applied to the minimum bias Monte 
Carlo, a rate of 5 • 10~^ candidates per interaction is observed, or approxi- 
mately 10"^ candidates per BX for the analyzed run. The rate of target single 
diffractive events (0.004 per BX, see Table 1) is also small compared to the ob- 
served rate. We conclude that the observed tracks cannot be due to hadronic 
interactions. 
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Xip(cm) X(cm) 



Fig. 5. a) The distribution of Xip for selected tracks (see text), b) The distribution 
of X-view impact points on the first VDS detector layer, for the same tracks used 
in part a. In both cases, the data is indicated by the histogram and the Monte Carlo 
by the points with error bars. 

Assuming that the observed tracks are indeed 5-rays, Equation 24 gives a lu- 
minosity estimate of 633 ± 28 ± 95 mb~^ per BX, to be compared with the 
luminosity estimate of 688 ± 35 mb~^ computed using the method described 
in Section 6. The agreement within errors lends further credence to the hy- 
pothesis that the tracks described in this section are 5-rays from the target 
and also serves as a cross-check of the method of Section 6. 



8 Summary and conclusions 



As previously noted the uncertainties affecting the total luminosity measure- 
ment are dominated by the systematic contribution, since each I A trigger run 
contains enough random trigger events to make the contribution from statis- 
tics negligible. In Table 6 we summarize the overall relative uncertainty on the 
total luminosity calculation (^^^). In the second column the uncertainty on 
Ka is given. This contribution depends on the present knowledge of the cross 
sections (see Table 1) and can in principle be improved in the future. The 
following three columns list individual contributions to the systematic uncer- 
tainty on detection which are combined according to Equation 15 to give the 
total detection uncertainties shown in the sixth column. 

When the method described in this paper is applied to the 2002 HERA-B 
minimum bias data taking period, the following integrated luminosities are 
obtained for each of the three target materials: 

C^^^ c = 405.8 ± 9.3 ± 17.9 fxh'^ 
£j„^^y. = 30.9 ±0.6 ±1.5 fxb-^ 
£^„i ^^ = 38.3 ±0.6 ±1.5 fib-^ 

where the first error corresponds to the uncertainty on Ka and second sum- 
marizes the remaining uncertainty mainly due to the HERA-B experimental 
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SKa 




(SC) 

\ C J Pni'i 


\ C J hkn 


(sct.A 


SCtot 


(SO 


c 


0.023 


0.039 


0.009 


0.018 


0.044 


0.050 


0.039 


Ti 


0.018 


0.042 


0.009 


0.023 


0.049 


0.052 


0.042 


W 


0.016 


0.032 


0.010 


0.019 


0.039 


0.042 


0.029 



Table 6 

Values of the contributions to the relative systematic uncertainty and overall relative 
systematic uncertainty on the total luminosity calculation ( ^^^°* ). The uncorrelated 




is given in the last column. 



conditions. The overall systematic uncertainty can then be obtained as the 
quadratic combinations of these two terms. 



The method for luminosity measurement described in this paper is based on 
the determination of the average number of interactions per bunch crossing, A, 
and on the knowledge of the total interaction cross section a (see Equation 2). 
The availability of a small fraction (few percent) of events acquired in parallel 
to the main stream of data with a completely unbiased trigger (pseudo-random 
trigger) was used to evaluate A on run by run basis. 

The measurement of A has been performed by exploiting the information from 
a variety of subdetectors, without the use of any dedicated device. This strat- 
egy allowed to perform consistency checks and to obtain a conservative deter- 
mination of the systematic uncertainties of the measurement. 



In case the three data sets will be combined to determine the A-dependence 
of a cross section, possible correlations between the systematic errors have to 
be taken into account. The correlated error is dominated by the uncertainty 
on \jnh and the background correction and is estimated to be ~ 3%. The 
correlation coefficients vary between 0.90 and 0.92 for pairs of wires. The 
uncorrelated part of the systematic uncertainty ( ^^^) can be found in the 

V '-tot ) uc 

last column of Table 6. 



The same strategy could be applied in future experiments such as those under 
construction at the LHC at CERN in Geneva, once the corresponding cross 
sections have been measured. 



^ At the time of publication of [23,24], the correlated error was estimated to be 
2%, rather than the updated and more accurate value (3%) presented here. A re- 
calculation of the values and errors reported in [23,24] using this updated estimate 
results in negligable changes. 
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